Metal-insulator transition in the one-dimensional Holstein model at half filling 
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We study the one-dimensional Holstein model with spin-1/2 electrons at half-filling. Ground state 
properties are calculated for long chains with great accuracy using the density matrix renormalization 
group method and extrapolated to the thermodynamic limit. We show that for small electron- 
phonon coupling or large phonon frequency, the insulating Peierls ground state predicted by mean- 
field theory is destroyed by quantum lattice fluctuations and that the system remains in a metallic 
phase with a non-degenerate ground state and power-law electronic and phononic correlations. When 
the electron-phonon coupling becomes large or the phonon frequency small, the system undergoes a 
transition to an insulating Peierls phase with a two-fold degenerate ground state, long-range charge- 
density-wave order, a dimerized lattice structure, and a gap in the electronic excitation spectrum. 
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A long time ago Peieridj suggested that a one- 
dimensional metal should exhibit an instability against 
a periodic lattice distortion of wave vector equal to twice 
the Fermi wave vector. Although this distortion in- 
creases the lattice elastic energy, it opens a gap in the 
electronic spectrum at the Fermi surface, lowering the 
electronic energy. Thus, the Peierls insulating state can 
be energetically favored over the metallic state. A wide 
range of quasi-one-dimensional materials, such as MX 
chains, charge-density-wave (CDW) compounds, conju- 
gated polymers and charge-transfer salts, have elec- 
tronic properties that are dominated or at least affected 
by the Peierls instability. These systems are often mod- 
eled by the one-dimensjonal Holstein model JL the Su- 
Schrieffer-Heeger modelfl or various spin-Peierlsa models. 

The Peierls instability is well understood in the static 
lattice (adiabatic) limit and within mean-field theory. 
An interesting and still controversial question is how 
the Peierls ground state is modified when quantum lat- 
tice fluctuations are taken into account. These quan- 
tum lattice fluctuations could have an important effect 
in most quasi-one-dimensional materials with a Peierls 
ground state because the lattice zero-point motion is of- 
ten comparable to the amplitude of the Peierls distor- 
tion.B Thus, this question has motivated several st\ 
of quantum lattice fluctuation effects in the-Hplsteir 
Su-Schrieffer-Heeger£3~t£l and spin-PeierlsE3 E3 models. 
In spinless fermion models and spin-Peierls models these 
studies have shown that the transition to a Peierls state 
occurs only when the electron-phonon coupling exceeds a 
finite critical value or when the phonon frequency drops 
below some finite threshold value. Thus, in these systems 
quantum lattice fluctuations destroy the Peierls instabil- 
ity for small electron-phonon coupling or large phonon 
frequency. In more realistic mQd£ls || w.ilb spin-1/2 elec- 
trons, however, previous studieaOoEjO'Ej have generally 
concluded that the ground state is a Peierls state for 
any finite electron-phonon coupling at finite phonon fre- 
quency, in qualitative agreement with mean-field theory. 

Here we consider the one-dimensional Holstein model 
with spin-1/2 electrons at half-filling. This model de- 



scribes electrons coupled to dispersionless phonons, rep- 
resented by local oscillators. It has as Hamiltonian 



(f) 



where and pi are the position and momentum opera- 
tors for a phonon mode at site i, c\ a and Ci yCr are creation 
and annihilation operators for an electron of spin a on 
site i, and m 
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The half-filled band case 



corresponds to a density of one electron per site. At first 
sight, there are four parameters in this model: the oscilla- 
tor mass M and spring constant K , the electron-phonon 
coupling constant a, and the electron hopping integral 
t. However, if phonon creation and annihilation opera- 
tors are denoted by b\ and bi, respectively, the Holstein 
Hamiltonian can be written (up to a constant term) 



H = cjJ24h-7Yl( b t + h) K-l) 

i i 

^ErHi/^ + ^/.+iJ , (2) 



where the phonon frequency is given by to 2 — K/M (we 
set Ti — 1) and a new electron-phonon constant is defined 
by 7 = a a with the range of zero-point phonon position 
fluctuations given by a 2 = ui/2K. We can set the param- 
eters t and a equal to 1 by redefining the overall energy 
scale and the units of phonon displacements. Thus, the 
properties of the Holstein Hamiltonian (2|) depends only 
on the two interaction parameters u> and 7. 

Mean-field theory predicts that the ground state of this 
model is a Peierls state for any non-zero electron-phonon 
coupling and w < 00. Early works based on strong- 
coupling perturbation theory and quantum Monte-Carlo 
simulations ,0 as well as variational calculations^ seemed 
to support this point of view. However, the quantum 
Monte-Carlo results were limited to small systems (up to 
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16 sites) and their interpretation relied on a questionable 
finite-size-scaling analysis. The strong-coupling pertur- 
bation theory is based on the formation of small bipo- 
larons in the j/lu — > oo limit, but it has been argued that, 
as the coupling 7 decreases, the bipolaron size becomes 
large and the strong-coupling picture breaks dowr£3. On 
the other hand, a functional integral calculation sug- 
gests that-the transition occurs at finite electron-phonon 
couplingt2l, but the accuracy of this approach is hard to 
estimate. Moreover, the static and dynamical properties 
of small clusters (up to six sites) show that there is a 
sharp crossover at a finite electron-phonon coupling from 
a quasi-free electron ground state to an ordered bipola- 
ronic ground state, which can be seen as a precursor to 
the Peierls ground state of the infinite system.Ej~Eil 

In this paper we discuss the ground state properties 
of the Holstein model of spin-1/2 electrons in the ther- 
modynamic limit. We demonstrate that quantum lat- 
tice fluctuations suppress the Peierls instability for small 
electron-phonon coupling or large phonon frequency. In 
this regime the ground state is unique, gapless and shows 
only power-law correlations between electron position 
and between phonon displacements. This ground state 
is similar to the ground state of the non-interacting sys- 
tem (7 = 0). When the electron-phonon coupling be- 
comes large or the phonon frequency becomes small the 
system undergoes a transition to an insulating Peierls 
phase, which is qualitatively described by mean-field the- 
ory In this regime the ground state is doubly degenerate, 
and there is a gap in the electronic spectrum, long-range 
CDW order and a dimerized lattice structure. 

Our results are based on density matrix renormaliza- 
tion group (DMRG) calculations SB DMRG is as accurate 
as exact diagonalization on small systems but can be ap- 
plied to much larger systems while maintaining very good 
precision. It has already been applied successfully to the 
study of the Peierls instability in quantum-J ^ttice s with 
spinless fermion or spin degrees of freedom .Ell 1 2 lr i There 
have not been any application of DMRG to models of 
spin-1/2 electrons coupled to phonons yet, because these 
systems are significantly harder to deal with due to the 
additional degrees of freedom and the larger amplitude of 
phonon displacements. For this work we have used an im- 
proved DMRG method for systems with boson degrees p£ 
freedom, which has been described in a previous work.LJ 
With this approach both the error due to the necessary 
truncation of the phonon Hilbert space and the DMRG 
truncation error can be kept negligible. The accuracy 
of this new DMRG technique has been demonstrated by 
comparison with many numerical and analytical methods 
for the polaron problem (a^ngle electron) in the one- 
dimensional Holstein model.13 l3 The maximum number 
of density matrix eigenstates m used in our calculations 
is 600, giving truncation errors from 10~ 7 to 10 -11 de- 
pending on the system size and parameters. The error 
in the ground state energy is estimated to be smaller 
than 10~ 5 t. The actual number of phonon states kept 
for each local oscillator ranges from 8 to 32 depending 



on the electron-phonon coupling strength. We have stud- 
ied open chains with an even number N of sites (up to 
100) and extrapolate results to the thermodynamic limit. 
Open boundary conditions are used because the DMRG 
method usually performs much better in this case than 
for periodic boundary conditions. 

In previous studies of the Peierls instability in the Hol- 
stein model the ground state symmetry was explicitly 
brokeg_aa. in the mean-field and adiabatic approxima- 
tions.Lrao Thus, the Peierls ground state was revealed 
by a lattice distortion (dimerization) 

= (-iy mp (3) 

and a CDW 

(m) = l + (-l) i m e (4) 

with m e ,m p 7^ 0, where (O) means the ground state 
expectation value of operator O, and m p and m e arc the 
phonon and electronic order parameter, respectively.!!] For 
m e , m p ^ 0, the ground state was two-fold degenerate 
[this degeneracy corresponds to the two possible phases of 
the oscillations (||) and (Q)]. Note that, as all eigenstates 
of the Holstein Hamiltonian satisfy 

(<ft> = J(W-l), (5) 
the order parameters are related by 

m p = ^m e . (6) 
K 

Our DMRG method gives an excellent approximation 
to the exact ground state of the Holstein model on a lat- 
tice of finite size. It is known exactly that the ground 
state of the half-filled Holstein model on a finite lattice 
is unique for lu ^= 0, implying that there is no degener- 
ate broken symmetry ground state at any finite electron- 
phonon coupling or non-zero phonon frequency.L^ In- 
stead, there is a quasi-degeneracy of the ground state 
when the-electron-phonon coupling exceeds a finite criti- 
cal valueEJ (this point will discussed in more detail later). 
Therefore, we always find (<&) = and (n.j) = 1 in our cal- 
culations. This property follows directly from the unique- 
ness of the ground state and the electron-hole symme- 
try, i.e., the invariance of the Hamiltonian (Q) under the 
transformation 

4 -> (-l) l <^ , q t -q t . (7) 

To observe the consequences of the Peierls instability we 
have to look at correlation functions. The most impor- 
tant ones for a Peierls state are the staggered charge den- 
sity correlation function 

C n (m) = (-l) m (( ni n i+m ) - 1) (8) 
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FIG. 1. Staggered charge density correlation function 
C„(rre) (solid line) and staggered phonon displacement cor- 
relation function C q (m) (dashed line) in the metallic phase 
for 7 = 0.4 (a) and in the Peierls phase for 7 = 1 (b). The 
distance m is calculated from the middle of an open chain of 
80 sites (a) and 40 sites (b), respectively. In both cases ui = 1. 

and the staggered phonon displacement correlation func- 
tion 



C q {m) = (-l) m ( mt+m ). 



(9) 



We have found that, for small electron-phonon coupling 
7 or large phonon frequency u>, both correlation functions 
decrease as a power-law m~P with 2 > j3 > as a func- 
tion of the distance m. An example is shown in Fig. |l](a). 
As the electron-phonon coupling increases or the phonon 
frequency decreases, the exponent (3 becomes smaller. 
For sufficiently large electron-phonon coupling or small 
phonon frequency the behavior of both correlation func- 
tions is completely different. As seen in Fig. [l](b), in 
this case both functions tend to finite values at large dis- 
tances, showing the existence of long-range order. 

It is not always possible to determine the presence or 
absence of long-range order in the thermodynamic limit 
from the correlation functions of a finite chain. A better 
approach is to compute the electronic and phononic static 
staggered susceptibilities defined as 



(10) 



0.06 



0.04 



0.02 - 



0.00 




0.00 0.01 0.02 0.03 0.04 0.05 

(N+1)- 1 




0.05 



FIG. 2. Electronic (circle) and phononic (square) staggered 
static susceptibilities as a function of the inverse chain length 
in the metallic phase for 7 = 0.4 (a) and in the Peierls phase 
for 7 = 1 (b). In both cases ui = 1. Solid lines are linear fits. 



and 



(ii) 



respectively. It is clear that both Xe and Xp vanish in 
the thermodynamic limit if there is no long-range order. 
For instance, both susceptibilities vanish as 1/N in the 
non-interacting limit (7 = 0). In Fig. ||(a) we show both 
Xe and Xp as a function of the inverse chain length for a 
weak electron-phonon coupling. Both quantities clearly 
tend to zero in an infinite chain. Thus, we conclude that 
there is no long-range CDW order nor lattice distortion 
in the ground state of the Holstein model for the param- 
eters (7 = 0.4,o; = 1) used in this example. On the other 
hand, it is clear that Xe and Xp remain finite for N — > 00 if 
there is long-range CDW order or a lattice dimerization, 
respectively. For instance, in the mean-field approxima- 
tion, one finds 



Xe 



Xp 



(12) 



Figure g(b) shows Xe and Xp as a function of the inverse 
system size for a relatively strong electron-phonon cou- 
pling. In this case, both susceptibilities remain finite for 
N — > 00 and thus, reveals the presence of a Peierls state 
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with long-range CDW order and lattice dimerization for 
the parameters considered in this example (7 = 1, u) = 1). 
Using Eqs. (^|) and (12), one sees that 



(13) 



in the mean-field approximation. It is possible to demon- 
strate that this relation holds for the exact ground state 
in several special cases, such as the adiabatic limit (to — > 
0) and the anti-adiabatic limit (lu — » 00). Although we 
can not prove the validity of (O) for the general case, 
our numerical results show that it is always satisfied 
(within numerical errors) in an infinite system. This sim- 
ply means that lattice dimerization and CDW are two in- 
separable features of the Peierls ground state. Therefore, 
we define a unique order parameter A as 



(14) 



where Xv an d Xe are the infinite system extrapolation 
of the ground state susceptibilities (|h]) and (|ll|) calcu- 
lated from DMRG simulations. If the ground state of 
the Holstein model is a Peierls state, one has A > 0, and 
otherwise A = 0. Obviously, this definition of A is just a 
generalization of the usual gap parameter of mean-field 
theory Amf, which is related to the other mean-field 
order parameters m e and m p by 



A M f = a\m p \ = -gWe 



(15) 



In the mean-field approximation the Peierls distortion 
opens a gap 2Amf in the electronic spectrum. It is some- 
times assumed that this relation between Peierls gap and 



n quantum lattice 
In such a case the 



order parameters remains valid w. 
fluctuations are taken into account, 
exact Peierls gap would simply be given by 2 A. How- 
ever, it is likely that the Peierls gap is more reduced by 
the quantum lattice fluctuations than the dimerization 
or CDW amplitudes and becomes smaller than the value 
2 A obtained from (|l4|). Unfortunately, calculating the 
optical gap of the Holstein model with a DMRG method 
is not possible yet.EJ To find how the appearance of the 
Peierls ground state correlates with a gap in the infinite 
system we have calculated the charge gaps 



E, 



2(^(1) -£o(0)) 



and 



E g2 =E (2)-Eo(P) 



(16) 



(17) 



where Eq(x) is the DMRG ground state energy with x 
electrons added to (x > 0) or removed from [x < 0) the 
half-filled band. In these definitions we implicitly use the 
electron-hole symmetry of the model at half filling, which 
implies that E (— x) = E (x). It should be noted that 
with these definitions the charge gaps incorporate lattice 
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FIG. 3. Gaps E g \ (circle) and E g 2 (square) vs. the inverse 
chain length in the metallic phase for 7 = 0.4 (a) and vs. 
the square of the inverse chain length in the Peierls phase for 
7 = 1 (b). In both cases u = 1. Solid lines are linear fits. 

relaxation effects occurring when the band filling is mod- 
ified. Therefore, E g \ and E g 2 are not always equal to 
the optical gap of the system. E g \ can be interpreted as 
the energy required to create a quasi-particle excitation 
made of an electron dressed by phonons. Similarly, 2E g 2 
represents the energy required to create a quasi-particle 
excitation which is a bound pair of electrons dressed by 
phonons, when such electron binding occurs (E g 2 < E g i). 
Otherwise, one expects E g2 ~ E g i . Figures ||(a) and (b) 
show both gaps for several system sizes. If there is no 
long-range order (A = 0) we find that the gaps extrapo- 
late to zero in the limit N — > 00 [Fig. 9(a)]. Therefore, we 
think that in this regime the system is still a metal, as in 
the non- interacting case (7 = 0). However, if the ground 
state of the infinite system is a Peierls state ( A > 0) , we 
find that both gaps extrapolate to a non-zero value in the 
thermodynamic limit [Fig. [|(b)]. For 7 = 1 and ui = 1, 
E g i = 0.82 and E g 2 = 0.18, which are much smaller than 
the value that one would anticipate from the amplitude 
of the Peierls distortion 2 A = 2.5. For comparison, the 
mean- field result for the same parameters is 2Amf = 3.1. 
This confirms that the quantum lattice fluctuations have 
a much stronger effect on the Pei«rls gap than on the am- 
plitude of the Peierls distortion^ Nevertheless, we have 
never found that either E g \ or E g 2 vanishes for N — > 00 
in the Peierls ground state. In small clusters, a sharp 
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FIG. 4. Lowest excitation energy ei (circle) and second 
lowest excitation energy £2 (square) as a function of the sys- 
tem size N in the metallic phase for 7 = 0.4 (a) and in the 
Peierls phase for 7 = 1 (b). In both cases u = 1. Solid lines 
are linear fits. 

drop of the Drude weight occurs simultaneously with 
the crossover to the ordered bipolaronic ground state.Eil 
Therefore, the opening of the electronic gap always seems 
to accompany the appearance of long-range order in the 
ground state and we conclude that a Peierls ground state 
is always an insulator. 

We have also analyzed the scaling of the lowest ex- 
citation energies e n — E n — Eq with the system size, 
where E n is the energy of the n-th lowest eigenstate of 
the Hamiltonian (|^) at half filling. In the phase without 
long-range order we have found that the e n decrease as 
a power-law for increasing system size and vanish in the 
thermodynamic limit, as seen in Fig. ^(a). These results 
confirm that in this case the infinite system has a unique 
ground state but is gapless; there is a continuous band of 
excitations starting from the ground state, as expected 
for a metal. In the Peierls phase, the energy difference 
£1 between the ground state and the first excited state 
is very small even in small chains and the other excited 
states have a much higher energy. Thus, the ground state 
appears almost degenerate in finite systems. Moreover, 
we observe completely different scalings for the e n . Fig- 
ure [|(b) shows that E\ decreases exponentially with in- 
creasing system size, while the energy differences between 
the two lowest eigenstates and the higher excited states 



FIG. 5. Order parameter A (circle) as a function of the 
electron-phonon coupling 7 for ui = 1. The solid line is the 
mean-field result Amf- 

remain finite in the thermodynamic limit. This shows 
that the ground state of the Peierls phase is two-fold 
degenerate in the thermodynamic limit. We have also 
checked that the order parameter A calculated for the 
first excited state tends to the same finite value as for 
the ground state in the thermodynamic limit. Therefore, 
both states are Peierls states with long-range CDW order 
and lattice dimerization, in qualitative agreement with 
mean- field predictions. The gap between the degenerate 
ground state and the other eigenstates also confirms the 
insulating nature of the system in the Peierls phase. 

Our results demonstrate that the ground state of the 
one-dimensional Holstein model for spin-1/2 electrons at 
half filling can be either a metallic state or an insulating 
Peierls state depending on the interaction parameters 7 
and lo. The system undergoes a quantum phase transi- 
tion between the metallic phase and the Peierls insulating 
phase at finite critical values 7 C and to c . In this aspect, 
the Holstein model for spin-1/2 electrons is similar to 
spin-Peierls and spinless fermion models. Unfortunately, 
DMRG simulations become less accurate and harder to 
carry out in the vicinity of the transition while, at the 
same time, the finite-size-scaling analysis requires more 
accurate results and larger system sizes. Therefore, de- 
termining the critical values j c and uj c for which this 
metal-insulator transition occurs demands a substantial 
amount of computer time and we have not attempted to 
draw a phase diagram. Nevertheless, we can show the 
evolution of the order parameter A as a function of the 
electron-phonon 7 for u> = 1 in Fig. [| We see that the 
transition to the Peierls state occurs around 7 = 0.8. 
This is in good agreement with calculations based on a 
functional integral approachJL9 which predicts 7 C ~ 1 for 
a slightly larger phonon frequency u = 1.1. As the adi- 
abatic and anti-adiabatic limits are usually investigated 
for finite values of the electron-phonon coupling constant 
A = a 2 /2K (= 7 2 /w with our choice of units), we show 
A as a function of the phonon frequency w for a fixed 
value A = 0.64 in Fig. One can see that our results 
converge to the exact adiabatic result for small u> and 
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FIG. 6. Order parameter A as a function of the phonon 
frequency uj for A = "f 2 /uj = 0.64. For uj — we show the 
exact adiabatic result. 

that the transition from the Peierls phase to the metallic 
phase occurs around w = 1. 

In summary, we have studied the ground state proper- 
ties of the one-dimensional Holstein model for spin- 1/2 
electrons at half filling using DMRG. We have shown 
that this system undergoes a transition from a metallic 
phase to an insulating Peierls phase at finite values of the 
electron-phonon coupling and of the phonon frequency. 
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